Decoherence and dissipation during a quantum XOR gate operation 
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The dynamics of a quantum XOR gate operation in a two-qubit system being coupled to a 
bath of quantum harmonic oscillators is investigated. Upon applying the numerical quasiadiabatic 
propagator path integral method, we obtain the numerically precise time-resolved evolution of this 
interacting two-qubit system in presence of time-dependent external fields without further approx- 
imations. We simulate the dissipative gate operation for characteristic experimental realizations of 
condensed matter qubits; namely, the flux and charge qubits realized in superconducting Josephson 
systems and qubits formed with semiconductor quantum dots. Moreover, we study systematically 
the quality of the XOR gate by determining the four characteristic gate quantifiers: fidelity, purity, 
the quantum degree, and the entanglement capability of the gate. Two different types of errors in 
the qubits have been modelled, i.e., bit-flip errors and phase errors. The dependence of the quality 
of the gate operation on the environmental temperature, on the friction strength stemming from the 
• system-bath interaction, and on the strength of the interqubit coupling is systematically explored: 

'^'n ' Our main finding is that the four gate quantifiers depend only weakly on temperature, but are rather 

\^ , sensitive to the friction strength. 
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iTi '. I- INTRODUCTION 

. The basic elements of quantum computation are logic quantum gates which represent manipulations of quantum 
' bits, |0) and |1), according to Boolean algebra. Any arbitrary complex logic operation can be build up of only a few 
basic gates {universal gates) [|l[ and one can show that almost every gate which operates on two or more qubits is a 
universal gate |^ . The explicit construction of quantum networks for elementary arithmetic operations then becomes 
Ch I possible upon appropriately combining such universal gates; see, for instance, Ref. |^ for the explicit construction of 
I i: the addition or the modular exponentiation. In turn, this permits the implementation of Shor's quantum factorizing 
■ O ■ algorithm Q in terms of elementary gates. Together with Deutsch's algorithm these two quantum algorithms are 
^ i: presently the most important examples which are known to be superior to their classical counterparts and which do 
justify the current efforts towards a technological realization of a quantum computer. 

In this work we concentrate on one such elementary gate, namely the quantum exclusive OR (XOR) gate. It 
is a unitary transformation which propagates an initial state l^'in) of a two-qubit system to a final state l^'out) = 
t^xORl^in)- Represented in the computational basis \bi) G {|00), |01), |10), |11)} {i = 1, ...,4), the XOR gate operation 
can be written as 
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UxoH = I 1 I ■ 
10, 

Since this operation inverts the state of the second qubit of the basis states if the first qubit is in the state |1), this 
operation is also called the quantum controlled NOT (CNOT) gate. The set of all onc-qubit gates together with the 
quantum XOR gate is universal, as has been demonstrated in Ref. 

The main impediment on the roadway to a working quantum computer is decoherence It disturbs the 

phase relation in a quantum superposition state and therefore is effective at the roots where the quantum computer 
is believed to have its most important advantage. Any realistic quantum computer will have some interaction with 
its environment which induces decoherence (decay of the off-diagonal elements of the reduced density matrix) and 
dissipation (change of populations of the reduced density matrix). Moreover, other sources for decoherence which are 
due to imperfect gate operations and to cross-talks of the qubits within a register need to be considered . 

Several previous works in the literature deal with the effect of decoherence in quantum information processing 
systems. Unruh |^ and Palma et al. |^ consider a model of a single qubit which is represented by the eigenstates of 
the quasi-spin operator cr^ and which couples to a bosonic environment via its (Tz-component. It describes appropriately 
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the dephasing (decoherence) but does not include population exchange (dissipation). Combining L non-interacting 
qubits of this type, they estimate the decoherence (in the limit of a large coherence length of the bath) to increase 
exponentially with the length L of the register. 

Dissipative effects (bit flip errors) are properly described by the so-called spin-boson model jl|-|l6|, where the qubit 
is represented by the o'a;-component of the spin-1/2, but the coupling to the bosonic bath is mediated by the a^- 
component of the spin-1/2 (note, that this refers to the localized representation). In this model, the bath also induces 
transitions between the two system eigenstates (bit flips) and - in addition to decoherence - energy is exchanged 
between system and bath. The general solution of the problem in terms of a generalized (non-Markovian) master 
equation for the entire reduced density matrix for an arbitrary initial preparation in presence of a static bias and also 
for a time-dependent driving has been given in Ref . [p^ . In our work, the assumption of a Markovian bath and of a 
weak system-bath interaction (Bloch-Rcdfield approach), which may restrict the validity of the master equation (see 
below), is not made. 

The previously discussed works concern the investigation of decoherence in single qubits or in a register of non- 
interacting qubits. Decoherence and dissipation in a system of interacting qubits has been studied only rarely. The 
dynamics of two coupled two-level systems has been investigated by Dube and Stamp ||l^ by means of a general 
model for coupled Josephson junctions, for coupled nanomagnets or for interacting Kondo impurities. Each two-level 
system is represented (in the tunneling representation) by the (Ta;-component of a spin-1/2. The two spins interact via 
their (Tz-components. Moreover, their (Tz-components couple to a bosonic bath. By use of real-time path integrals the 
dynamics of the relaxation process is determined. Although no specific problem of quantum information processing 
is investigated, this is the first work where two interacting spins in a dissipative bath have been considered. 

A similar model has been studied by Governale, Grifoni and Schon [|8|. Two biased spin-1/2 systems interact via 
their (jy-components which is the appropriate coupling for Josephson junction charge qubits (see below). Moreover, 
their (Tz-components couple either to the same or to different bosonic baths. Applying the widely used Bloch-Redfield 
formalism, the time evolution of the populations of the logical states is evaluated. This model describes dissipation 
being caused by fluctuations in voltage sources in Josephson junction charge qubits (see below). However, no specific 
quantum information operation has been considered. 

A two-qubit quantum gate for quantum information processing in coupled quantum dots has been investigated in 
Refs. ||l^,|l^. Two spin-1/2 systems are coupled using a time-dependent Heisenberg-type interaction. Moreover, a 
coupling of the spins to a bosonic bath has been taken into account. By solving the quantum Liouville equation in the 
limit of weak system-bath coupling (Born-Markov approximation) for the reduced density operator, the purity and 
the fidehty of the swap operation J7swap|ij) = — 0, 1) is calculated as a function of time. However, the authors 

consider the time-evolution of the quantum system after the swap operation has been completed. The same is true for 
the XOR gate operation in Ref. ||ll[], where, additionally, a further assumption has been made: The pulse sequence 
to realize the quantum XOR consists of four pulses of the external fields. Each pulse is taken to be constant over the 
corresponding time interval. To obtain the solution over the entire time span within the Born-Markov approximation, 
it is necessary to assume a finite time interval between the single pulses. This is required because the Born-Markov 
approximation is known to violate positivity of the reduced density operator at short transient times ||2C| , |2l| . This 
additional time span (pulse-to-pulse time) has been taken as three times the switching time interval. This leads to 
an extension of the computation time which is only due to formal mathematical reasons and which deteriorates the 
quality of the gate operation. Moreover, a systematic study of the dependence of the gate quantifiers on the relevant 
parameters has not been given. 

In this work, we investigate systematically the XOR quantum gate in presence of an interaction of the qubits 
with their environment. Thereby, we take into account the full time-dependence of the external fields which induce 
the XOR operation without invoking further approximations on the system Hamiltonian. In particular, we use the 
numerical ab-initio technique of the quasiadiabtic propagator path integral (QUAPI) p^ ] (for other applications, see 
also Refs. [^,Q). This numerically precise iterative real-time path integral method does not suffer from the above 
mentioned problem of lacking positivity. In order to realize the logic XOR operation in physical systems, we introduce 
a generic model Hamiltonian which is suitable for studying the XOR operation on a very general and idealized level. 
We determine the quality of the gate by calculating the four characteristic gate quantifiers introduced by Poyatos, 
Cirac and ZoUer P5[ : namely, the (i) purity, (ii) fidelity, (iii) quantum degree, and (iv) entanglement capability. To 
that end, we consider two important types of computational errors, i.e., phase errors and bit fiip errors. The former 
can be modelled by coupling the CTz-component of each spin to the bath while the later is induced by coupling the 
CTaj-component of each spin to the bath. We are mainly interested in the quality of the gate operation during its 
time evolution and, most importantly, right after it has been completed. Our analysis covers a wide field of physical 
systems and we choose parameter sets which mimic the realistic physical situation for flux qubits or charge qubits 
proposed in superconducting hybrid systems as well as charge and spin states in coupled semiconductor quantum dots 
(see below). 

So far, we have discussed theoretical aspects of quantum information processing. However, those refined and highly 



2 



elaborate concepts face the question of how they can be implemented in experimental hardware. Several proposals to 
build a quantum information processor exist. Prominent candidates are, for instance, atoms in optical cavities, ions in 
linear or Paul traps interacting with laser beams, or nuclear spins in an NMR liquid |^. Although the experimental 
techniques in those fields of research are currently most advanved, the problem of upscaling of a quantum computer can 
seemingly only be solved within condensed matter systems which can be embcded in an electronic circuit. Promiment 
systems for condensed matter qubits are flux states of a superconducting quantum interference device (flux qubits) 
prf (see also |[T^]), charge states of superconducting islands with Josephson junctions [charge qubits) p^ , |2^ , and 
spin ^1 29| or charge |30|] states in ultrasmall coupled semiconductor quantum dots (quantum dot qubits). Moreover, 



several realizations of qubits in nuclear 32 1 and electronic |3^,|3^ spins in semiconductor nanostructures have been 
proposed. In this work, we concentrate on typical experimental situations for flux and charge qubits in Josephson 
junction devices and for charge qubits in coupled quantum dots. 

The paper is organized as follows: In the subsequent Section ||, we introduce a generic model as a starting point 
for the quantum XOR operation including the interaction with the environment. In Section |l|, we present a brief 
review on the numerical technique of the quasiadiabatic propagator path integral (QUAPI) which we employ in the 
following. In order to determine the quality of the decoherent XOR gate, we use four quantifiers which are introduced 
in Section The results and the conclusions are presented in Sec. ^ and Sec. VI, respectively. 

II. A GENERIC MODEL FOR THE QUANTUM XOR GATE 

A. The coherent XOR operation 

The quantum XOR gate is a two-qubit operation which can be modelled by two coupled spin-1/2 systems represented 
by the Pauli operators aj — (aj,aj,aj)'^,j = 1,2. The two logical states of each qubit are represented by the two 
eigenstates of the cr^-component of each spin, i.e., \0)j = \ 1)j and = \ We assume that the single qubit as 
well as the coupling between the two qubits can be controlled by switching on (local) external fields, for instance, 
magnetic fields. This system can generically be described by the Hamiltonian 

h ^ 

Hxonit) = - 2 ^ + "t'^k ■ (2) 

where — (ctJ + ia^)/2. Moreover, Bj(t) = (Bj (t),0, Bj (t))'^ , j — 1,2 arc time-dependent coupling strengths 
(with the dimension of a frequency) arising from local time-dependent external fields at the site of the spin j in 
longitudinal (z-) or transverse (x-) direction. In Eq. (|2|), the coupling between the two qubits is assumed to be 
symmetric; furthermore it should be controllable from the outside leading to a time-dependent interaction strength 
J(t). The particular form of the interaction in Eq. (^) is only one example. We note that this generic model does not 
account for the particular details of a physical realization of qubits in a specific condensed matter system. For each 
individual system, such as flux qubits or charge qubits, the Hamiltonian looks different in detail. In particular, the 
coupling term between the two qubits takes different forms. However, the differences to our generic model in Eq. (H) 
are of minor influence only. The general physical behavior will be similar such that our generic model serves as an 
archetype. 

The quantum XOR gate (pi) can be obtained by a sequence of one- and two-qubit operations according to [p2| 



^xo« = aj(|)af(-|).,J(-.,t,.(-|)af(-|) 



x..(f)</r(-f)^i(-|), (3) 



where 



TJX/Z 

j 



(a) = exp (i^c^J^"") , j = 1, 2 , 
UMP) = exp {iP(ata^ + a^a+)) (4) 



are the propagators over the single time-intervals with the external fields in the Hamiltonian, Eq. (^), switched on and 
off in the following way: In order to attain this propagator, a pulse sequence of the external fields is necessary. For 
simplicity we assume throughout this work, that the pulses are switched instantaneously on and off and are constant 
over the time span toff — ^on they are on. This induces time-dependent interaction strengths B(t) = B[Q(t — ion) — 
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0(i — ioff)] with B = Bj^', J = const, and with Q{t) being the Heaviside function. Furthermore, we assume that 

both spins are equal and experience local fields of equal strength. This implies Bi^^ ~ Bj^^ = B^/^. The angles a 
and P in Eq. (|^) are related to the actual physical propagation time t according to 

a = B^/^t und (3 = Jt. (5) 

The switching times then follow as ti = 7r/(2B'^),i2 = h + 7r/(2J),t3 = i2 + 7r/(2B^),<4 = ts + 77/(2,7), = 
ti + n/{B^), = ts + it/{2B^) and txoR = ^6 + ''■/(2i?^), where ixOR denotes the total time elapsed during the full 
XOR gate operation. An example of this pulse sequence is sketched in Fig. |^ for the case of B'^ = B^ = J. The 
coupling constants are given in units of B^ while the time is scaled in units of {B^')~^. One immediately observes 
that the computation time txOR is extended if the coupling energies are decreased. We note that the assumption of 
rectangular pulses is not required by the numerical technique we use and is made here only for the sake of simplicity. 
We could also consider other shapes of the pulses which are more realistic for specific physical systems, and especially, 
we could consider imperfect switching processes as well; the latter would constitute a further source of decoherence. 



B. Interaction with the environment 



We model the interaction of the qubit system with the fluctuating environment by a Hamiltonian, in which i?xoR(0 
is coupled to a bath of harmonic oscillators, i.e.. 

Hit) - i?xoR(i) +Hb+ H^I: , (6) 

with 

^ / i\ 
HB=Y,^Aa^(^i + ^] . (7) 

Here, a^{aj) denotes the creation (annihilation) operator of the j-th bath oscillator with frequency uij. Since we 
want to investigate the role of bit-flip errors as well as phase errors we include in our model two different types 
of interactions: On the one hand, the cr^-components of the spins couple to the fluctuating environment and the 
populations of the qubit states are disturbed (bit-flip errors). On the other hand, phase errors are generated by a 
coupling of the cr^-components of the spins to the environmental noise. This is conveniently modelled by the form 

TjX/z / X/Z x/z\ Xjz, + N 

Hint =2^1 +"^2 jz^«/ (aj + aj), (8) 

where k^^^ denotes the coupling strength of the j-th oscillator to the system and where the superscript {x/ z) denother 
either the one or the other kind of interaction. We note that we assume here, to keep things simple, a coupling of the 
two spins to the same bath. This implies that the spins are effectively coupled to each other via the bath. A coupling 
of the spins to different (uncorrelated) baths could be readily incorporated in the numerical QUAPI technique (see 
below) . 

To study the dynamics of this system, we have to specify the initial conditions. Throughout this work, we assume 
that the density operator W{t) of the entire system-plus-bath at initial time t — factorizes according to 

W^(0) = ps(0)®pB. (9) 

PsiO) is the density operator of the system at time t = and pB = -^b ^ cxp {—H-q/ [k-QT)) is the canonical equilibrium 
distribution of the (decoupled) bath at temperature T. Moreover, Z-q = trexp(— i/B/(fcB2^)) and k-Q denotes the 
Boltzmann constant. 

The influence of the bath is fully characterized [Q by the spectral density 

N 

V^/^uj) = 2^Y.^nf)^5{u~u,) (10) 

which assumes a continuous form if the number N of oscillators approaches infinity. Throughout this work, we apply 
an Ohmic spectral density with an exponential cut-off, i.e., 

r-/^(c^) = 7-/^0. exp(-c^K) , (11) 
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where the dimensionless damping parameter 7^/^ characterizes the strength of the interaction with the environment. 
This spectrum mimics the environmentally induced fluctuations in the external circuit which supplies flux through the 
SQUID loops in the flux qubits • Moreover, background charge fluctuations in the voltage sources in Josephson 

charge qubits fl^ , ^^ also lead to an Ohmic impedance R. Similarly, electronic states in coupled quantum dot qubits 
experience an Ohmic environment, either for the spin 0,M or for the charge [Brt degrees of freedom. 



III. NUMERICAL AB-INITIO TECHNIQUE: QUAPI 

In order to describe the dynamics of the two-qubit system of interest it is sufficient to consider the time evolution 
of the reduced density operator 

-1/ 



p{t) ^tvi,,thU{t,0)W{0)U-\t,0) , 
U{t,0) ^Texpl^--i/h H{t')dt' 



(12) 

Here, U {t, 0) is the propagator of the full system-plus-bath and T denotes the time-ordering operator. Moreover, 
trbath means the partial trace over the harmonic bath oscillators. Due to our assumption that the bath is initially at 
thermal equilibrium and decoupled from the system, see Eq. (^), the partial trace over the bath can be performed. 
We denote the matrix elements of the reduced density matrix in the computational basis with pij{t) = {bi\p{t)\bj) 
and rewrite them according to Feynman and Vernon pit] as 



with the propagator G given by 



PijW ^ X! G'ij,™„(t, 0)p„„(0) , (13) 



G^J,mn{t,0)= J VxVx' A[x]A* [x']J'fv[x , x'] . (14) 



The functional A[x] denotes the probability amplitude for the free system to follow the path x{t) and J^pv[x, x'] 
denotes the Feynman- Vernon influence functional (see Ref. for details). The functional integrations in Eq. 
(H) extend over paths with endpoints x(0) = x^, x{t) = Xi, x'{0) — Xn and x' [t) ~ Xj which belong to the initial and 
final states, PmniO) and Pij{t), respectively. 

The technique which we use to calculate the reduced density operator Eq. (^3|) is the iterative tensor multiplication 
scheme derived for the so-called quasiadiabatic propagator path integral (QUAPI). This numerical algorithm was 
developed by Makri and Makarov |^ within the context of chemical physics. Since its first applications it has been 
succcsfully tested and adopted to various problems of open quantum systems, with and without external driving 
p2|-p4t. Because the details of this algorithm has been extensively discussed previously in the literature p2|-p^, we 
only mention those prominent ingredients which arc of importance for our work: 

(i) Symmetric splitting of the short-time propagator: To obtain a numerical iteration scheme, we discrctize the 
time interval [0, t] into J\f steps At, such that tk = kAt and split symmetrically the full propagator over one time step 
U{tk+i,tk) in Eq. ( p^ ) according to the Trotter formula into a system and an environmental part: 

U{tk+utk) « eM~^HBAt/2h)Us{tk+utk)cxp{-iHBAt/2h) , 



k + l 



Usitk+utk)^Texp^~-J dt'HxoR{t')j . (15) 

The neglect of higher order terms of the propagator in Eq. ([l5| ) causes an error of the order of At^. The short-time 
propagator Us of the bare system is given by the corresponding exact system propagators in Eq. (^) over a time step 
At. At this point, we emphasize that this method is not plagued by the problem of lacking positivity of the density 
operator at short times, as it is the case for the usually employed master equation approach in the Born-Markov limit 
| pT| , p^ . The exact coherent dynamics of the bare system enters and the decomposition of the short-time propagator 
according to Eq. ( p^ ) is valid for any arbitrary short time. 

(ii) The interaction with the bath induces correlations among the paths (memory) which are described by the 
influence functional in Eq. ([l^). As long as the temperature of the Ohmic bath is finite, these correlations decay 
exponentially fast with increasing time [H. This motivates to neglect such long-time correlations and to break up 
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the influence kernels into smaller pieces of length K At, where K denotes the number of time steps over which the 
memory is fully taken into account. 

The two strategies in (i) and (ii) are countercurrent: In step (i) a small time step At is desirable in order to 
minimizes the error due to the neglected higher order terms in the propagator. On the other hand, in (ii) a large time 
step is wanted in order to take a long memory range into account. A compromise between those two errors has to be 
found in practice by applying the principle of minimal sensitivity to adjust the two parameters At and K , see 
discussion below. 

(iii) The third important ingredient is the appropriate choice of basis representation of the problem. For the 
algorithm it is required to iterate the dynamics in the eigenbasis of that system operator which couples to the bath. 
Then the influence functional in Eq. (^4|) can be evaluated in terms of the eigenvalues of the coupling operator. In 
problems where the coordinate of a quantum particle in a continuous potential is damped, the continuous position 
operator turns into a discrete set of position eigenvalues. Hence, this representation has been termed discrete variable 
representation (DVR). 

Bit- flip errors: If the (Tj;-components of each spin couple to the bath, see Eq. (^), the eigenbasis of the coupling 
operator is determined by {ai\ {a'f + erf) /2|aj) = 5ij with Ai = 0, A2 = —1, A3 = 1 and A4 = 0. A basis rotation of 
the computational basis with basis states \bj) has to be performed according to 

4 

(16) 

with the transformation matrix 




(17) 



Phase errors: For the second case that the cr^-components of each spin couple to the bath in Eq. (|^), the system part 
of Hamiltonian is already diagonal in the computational basis, i.e., {hi\ {af +(^2) /-^l^j) = \ with Ai = 1,A2 = 
0, A3 = and A4 = —1. No additional basis transformation is necessary. 



IV. CHARACTERISTIC GATE QUANTIFIERS 



In order to quantify the quality of the quantum gate, we use four global parameters which have been defined by 
Poyatos, Cirac and ZoUer (i) the gate fidelity JF, (ii) the gate purity 7^,(111) the quantum degree Q of the gate, 
and (iv) the entanglement capability C of the gate. These four quantifiers can be calculated once the reduced density 
operator p in Eq. ( p^ ) is determined. To this end, 16 unentangled input states 1^';^):^ = 1, 16 are defined according 
to \^a)Mb)2 (a,6=l,...,4), with 1^1) = |0),|V'2) = = (|0) + |l))/V2and 1^4) = (|0) + z|l))/\/2 . They form 

one possible basis set and span the Hilbert space for the superoperator VxoR where p(txoR) = V'xorp(O), see Ref. 
p5[ for details. Moreover, these basis states are chosen to be unentangled states in order to avoid the application of 
a preceding two-qubit gate for the preparation of the system state. 

The gate fidelity T is defined as the overlapp between the propagation with the ideal propagator ?7xoRi Eq. (|^), 
averaged over all 16 initial states according to 

1 

^ = E(*in|f^XORPxORt^XOR|*L> (18) 

with p^Qj^ = p(txoR) with initial condition p{0) = |*in)(*inl- 
In a similar way, the purity V is defined as 

16 

^=I^Et'^KoR)'- (19) 
This quantity is proportional to the (negative) linearized entropy and reflects the effects of decoherence. 



6 



The third quantity, the quantum degree Q of the gate, is defined as the maximum of the overlap of all possible 
output states stemming from unentangled states and of all maximally entangled (Bell) states \'^^^),k = 1, ...,4. In 
formal terms, this implies 

S = max«Jp^oi,K,). (20) 

The purpose of this parameter is to quantify the notion of nonlocality. Bennett and co-workers | ]35[ | have shown 
that all those density operators which have an overlap with a maximally entangled state being larger than the value 
(2 + 3\/2)/8 « 0.78 are non-local, i.e., violate the Clauser-Horne-Shimony-Holt inequality p5[ . 
Obviously, T — \,V — 1,Q = 1 denote the ideal gate operation. 

The fourth quantifier is the entanglement capability C of the gate. It denotes the smallest eigenvalue of the 
partial transposed density matrix |36| which is determined from PxoR unentangled input states PxoR 

characterizes an entangled state if and only if the smallest eigenvalue of the partial transposed density operator is 
negative. The ideal operation has an entanglement capability of C = —0.5. 



V. RESULTS 



Having determined the reduced matrix in Eq. (|ij) by the iterative QUAPI algorithm, we investigate the influence 
of the interaction with the environment systematically. Therefore, we assume that the two qubits are identical and 
experience external fields of the same strength, i.e., B'f — _Bf — and El = i?| = B^. Moreover, we introduce 
the following dimensionless parameters: We scale the quantities with respect to the characteristic energy scale of the 
single qubit which is given by the energy splitting hB^ of the single qubit. This in turn defines a time scale {B^)~^. 
Consequently, the temperature is given in units of hB^ /k^ (Note that 7^/^ is already dimensionless). For all following 
results, we have used a cut-off frequency of Wc = 505^ in (O). 



A. Time resolved quantum XOR operation 



We first illustrate the time-resolved dynamics of a generic XOR operation. To this end, we determine the populations 
of the four states of the computational basis as a function of time, i.e., Pij{t) := {ij\p{t)\ij) with i, j = 0, 1 for the 
initial condition p(0) — |11)(11|. We choose the pulse sequence sketched in Fig. |l|with parameters set to B^ ~ B^ and 
J = B^. Moreover, we choose for illustrative purpose a rather high temperature of T = O.lhB^ /ku- Fig. || depicts 
the results for the three different cases of (i) no damping, 7^/^ — 0, (solid line), (ii) bit flip errors with 7^ = 0.01 
(long dashed line), and (iii) phase errors with 7^ = 0.01 (dotted line). The switching times tj are equal to multiples 
of 7r/2 for this special case of equal energies. 

The iterative QUAPI algorithm possesses two free parameters which have to be properly adjusted. We fix the 
number K of memory time steps and the length At of each time step according to the principle of minimal sensitivity 
p^ . By applying this method, we obtain the values At — 0.15(i?^)~"'^ with K — 2 (not shown). 

As one observes, the final state of the ideal operation (7^/^ = 0) is j^I') = |10). The deviation of the dynamics in 
presence of decoherence and dissipation from the ideal case is clearly visible. 



B. Quality of realistic quantum XOR operations 



In order to study realis tic ph ysical situations, we ap ply ou r method to three different condensed matter systems, 

charge qubits (set II) 



namely flux qubits (set I) |27 Ij], charge qubits (set II) [12[|2^], and qubits realized in coupled semiconductor quantum 
dots (set III) |l^,^^. We choose typical parameter sets published in the literature. They are summarized in Table |. 
The QUAPI parameters are again determined by the principle of minimal sensitivity for K = 3. We obtain for the 
(bit flip errors) for set I (in units of (B^)-i): At = 0.06, for set II: At = 0.2, and for set III: At = 0.02, 



case of H?^^ 



and for the case of H[^^ (phase errors) for set I: At = 0.013, for set II: At = 0.08, and for set III: At = 0.01. 



1. Dependence on temperature 



The dependence of the four characteristic gate quantifiers on the bath temperature T is depicted with Fig. ^ In 
panel a.) the influence of the random bit flips are investigated while the panel b.) depicts the effect of phase errors. 
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The damping constant for the bit flip errors is chosen to be 7^ = 10 ® and for the phase errors we seet 7^ = 10 

First, one observes that all results depend only weakly on temperature. Extrapolating the results to zero temper- 
ature indicates the influence of the zero point oscillations of the harmonic bath oscillators. Second, the results for 
!F,V and Q do in no case exceed a value of 0.999 85 for bit flip errors and 0.975 for phase errors. This demonstrates 
that even smaller strengths of the coupling to the environment than 7'^ = 10^^ or 7^ = 10^^ are necessary in order 
to obtain a desired value of 0.999 99 0! 

Furthermore, we observe that the flux qubits (set 1) are the least sensitive to bit flip errors while charge qubits 
(set II) and quantum dot qubits (set III) both are performing worse, see Fig. || a.). On the other hand, quantum dot 
qubits (set III) are most immune from phase errors, see Fig. ^ b.) while the superconducting systems assume worse 
results. 



2. Dependence on the damping strength 

As shown in the preceding section, the quality of the gate operation cannot be improved by lowering the temperature 
of the environment. On the contrary, the shielding of the qubit system against external noise is most important. This 
is demonstrated by the results for the dependence of the four gate quantiflers on the coupling constants 7^/^ depicted 
in Fig. |. 

We find again that the results for the flux qubits (set I) approach the ideal values olJ-=V = Q= \ and C = —0.5 
first upon decreasing the friction strength, see Fig. |j a.). However, the flux qubits are more sensitive to phase errors 
from which the quantum dot qubits (set III) are seemingly most immune, see Fig. ^ b.). Note that in Fig. ^the lower 
bound of Q ~ 0.78 for the Clauser-Horne-Shimony-Holt inequality is indicated by the horizontal dotted-dashed line. 

For all three systems we find that the damping strength has to be less than 7^/^ « 10~^ in order to serve as suitable 
candidates for a quantum information processor. 

3. Dependence on interqubit coupling strength 

In the remaining part we address the dependence of the quality of the gate operation on the strength J of the 
interqubit coupling. Physically, one expects that a coupling strength which is comparable to the characteristic qubit 
energy (J « i?^) would yield the best results. This is confirmed in Fig. ^jfor the two different sets I and III (set II is 
equivalent to set I). The panel a.) illustrates the results for the bit flip errors with 7^ = 10~^, and correspondingly 
panel b.) for the phase errors with 7^ — 10^*. As one can observe, the superconducting qubits (set I) can perform 
with a smaller interqubit coupling than the quantum dot qubits (set III). Moreover, the use of an interqubit coupling 
larger than J w 10^^ is advantageous. The horizontal dotted-dashed fine in Fig. ||b.) indicates the lower bound 
of Q w 0.78 for nonlocality. 

VI. CONCLUSIONS 

In this work we have shown that the numerical quasiadiabatic propagator path integral method (QUAPI) of Makri 
and Makarov provides an appropriate method to investigate decohering quantum information processes which involve 
time-dependent Hamiltonians in presence of a coupling to an external environment. We have applied this iterative 
algorithm to the example of the quantum XOR gate operation and have obtained the full time-resolved evolution of 
the two-qubit system in presence of time-dependent external fields. No further approximations on the time evolution 
of the gate operation such as a Markovian evolution or extended time spans of the gate operation have been invoked. 

We have investigated the quality of the quantum XOR operation by numerically determining four characteristic gate 
quantifiers, i.e., the fidelity J^, the purity P, the quantum degree Q, and the entanglement capability C of the gate. We 
have simulated the gate operation for three parameter sets which mimic realistic experimental situations for condensed 
matter qubits, such as flux and charge qubits in superconducting Josephson systems and qubits in semiconductor 
quantum dots. Thereby, we have suceeded to investigate systematically the quality of the gate operation as a function 
of the environmental temperature T, the friction strength 7, and the strength J of the interqubit coupling. Two 
different types of errors in the qubits have been modelled: (i) bit-flip errors and (ii) phase errors. We have elucidated 
how the different physical setups perform under these conditions. As major findings we establish that the quality 
of the gate depends only weakly on temperature but rather strongly on the friction strength. Moreover, we have 
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illustrated that the interqubit coupling strength plays an important role and should not be smaller than O.IEq with 
-Bo being the typical energy scale of the qubit system. 

In order for a quantum information processor to operate optimally the decoherent influence of the environment 
needs to be suppressed. Therefore, three different approaches are currently discussed: These are the techniques of 
quantum error correction, fault tolerant quantum computation, and entanglement purification pGf . The general idea 
common to all three methods is to use for quantum information processing only a small subset of a larger set of 
entangled ancilla qubits. Although these ideas are very promising for small register lengths, the techniques become 
increasingly difficult if one attempts to realize large qubit registers in physical systems. Moreover, one has to keep 
track of the quantum state of the environment. Whilst these requirements are seemingly feasible for quantum optical 
information processing systems [^7| , they appear insufficient for condensed matter systems with their characteristic 
huge number of environmental degrees of freedom. 

An alternative approach consists in minimizing the occurrence of errors by controlling decoherence via the ap- 
plication of tailored time-dependent external fields to qubit systems p5| , |3^ -p[. For example, the application of a 
time-dependent periodic external field can induce a Floquet spectrum with degenerate quasienergy states JT5| , p8| , ^ 
or it can move the qubit out off resonance with certain bath modes thereby reducing decoherence. The suitability 
of such schemes to a quantum gate operation, however, remains to be demonstrated. 
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Set 






J 


T 




J/B'' 


T/B' 


I: Flux qubits 


0.5 K 


50 mK 


25 mK 


25 mK 


0.1 


0.05 


0.05 


II: Chaxge qubits 


1 K 


100 mK 


5 mK 


50 mK 


0.1 


0.005 


0.05 


III: Quantum dot qubits 


1 meV 


1 meV 


0.05 meV 


125 mK 


1 


0.05 


0.01 



TABLE I. Parameter sets which have been used for the simulations of the quantum XOR gate. They mimic typical experi- 
mental situations in three important solid-state qubit systems, i.e., flux qubits and charge qubits in superconducting Josephson 
devices (Set I and II) and spin and charge qubits in ultrasmall semiconductor quantum dots (Set III). 
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Bi'(t) I 1 

Bz'Ct) I I I 1 

Bi\t) I 1 

Ba'Ct) I 1 I 1 

J(t) I 1 I 1 

tj t2 t3 t4 txoR 



FIG. 1. Schematic view of a pulse sequence necessary to generate the quantum XOR gate. The parameters are set to 
= = J = const.. The frequencies are given in units of while the time is scaled in units of {B^)~^ . The switching 
times tj are given in the text and are in this case equal to multiples of 7r/2. 
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XOR gate: T=0.1 , to^=50 
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FIG. 2. Time resolved dynamics of the quantum XOR operation for the case of equal energies, i.e., = = J. Depicted 
are the populations Pij{t) — {ij\p{t)\ij) as a function of time for the initial condition p(0) = |11)(11| for three different cases of 
(i) no damping, 7^/^ — 0, (solid line), (ii) bit flip errors with 7^ = 0.01 (long dashed line), and (iii) phase errors with 7^ = 0.01 



(dotted line). The time is scaled in units of (B^ 
frequency to Wc = 505^. 



Moreover, we set the temperature to T = O.lhB^ /kB and the cut-off 
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Dependence on temperature 

a.) Bit flip errors y'=W^ 
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b.) Phase errors y^=10 
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FIG. 3. Dependence of the fidelty J^, the purity V, the quantum degree Q, see - |2^) and the entanglement capabihty C 
on temperature T for bit flip errors (Fig. ^ a.) and phase errors (Fig. ^ b.). The temperature is scaled in units of HB^/kB- 
The qubit parameters are given in Table |. The damping constant for the bit flip errors is set at 7^ — 10~^ and for the phase 
errors at 7^ — 10""*. 
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Dependence on friction strength 

a.) Bit flip errors 
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FIG. 4. Dependence of the four gate quantifiers on the dimensionless friction strength. Upper panel a.): Bit flip errors (7^), 
lower panel b.): Phase errors (7^). The lower bound of Q « 0.78 for the Clauser-Horne-Shimony-Holt inequality is indicated 
by the horizontal dotted-dashed line (see text) . For the remaining parameters, see Table |. 
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Dependence on interqubit coupling 

a.) Bit flip errors J'=^0'^ 
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b.) Phase errors y^=10 




FIG. 5. Dependence of the gate quantifiers on the strength J of the interqubit couphng. Depicted are the resuhs for 
parameter set I and III from Table |. The upper panel a.) shows the results for the bit-flip error with 7^ = 10~^ while the 
lower panel b.) depicts the results for the phase error with 7^ — 10~*. The horizontal dotted-dashed line marks the lower 
bound of Q ~ 0.78 for the Clauser-Horne-Shimony-Holt inequality. 
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